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Abstract 

The infrared absorption spectrum of the protonated water dimer (HsO^) is simulated in full 
dimensionality (15D) in the spectral range 0-4000 cm^^. The calculations are performed using the 
Multiconfiguration Time-Dependent Hartree (MCTDH) method for propagation of wavepackets. 
All the fundamentals and several overtones of the vibrational motion are computed. The spectrum 
of H502^ is shaped to a large extent by couplings of the proton-transfer motion to large amplitude 
fluxional motions of the water molecules, water bending and water-water stretch motions. These 
couplings are identified and discussed, and the corresponding spectral lines assigned. 

The large couplings featured by HsO^ do not hinder, however, to describe the coupled vibra- 
tional motion by well defined simple types of vibration (stretching, bending, etc.) based on well 
defined modes of vibration, in terms of which the spectral lines are assigned. Comparison of our 
results to recent experiments and calculations on the system is given. The reported MCTDH 
IR-spectrum is in very good agreement to the recently measured spectrum by Hammer et al. 
[JCP, 122, 244301, (2005)]. 



* e-mail: |oriol.vendrell@pci.uiii-heidelberg.de| 
^e-mail: gatti@univ-niontp2.fr 



e-mail: Hans-Dieter.Meyer@pci.uni-heidelberg.de 



I. INTRODUCTION 



The understanding of the hydrated proton in bulk water and water-containing systems 
is of importance to many areas in chemistry and biology. Much effort has been directed 
in recent years towards obtaining a better understanding of the excess proton in water and 
important contributions have appeared due to the advance in experimental and theoretical 
techniques [l, 2, sl, JQ- 

The importance of the hydrated proton and the amount of work devoted to the problem 
contrast with the fact that the smallest system in which a proton is shared between water 
molecules, the protonated water dimer (HsO^), is not yet completely understood, and an 
explanation of some important spectral signatures and an overall picture of the dynamics 
of the cation is lacking. Accurate infrared (IR) spectroscopy of protonated water clusters 
prepared in the gas phase has become possible in recent years , opening the 

door to a deeper understanding of these systems and the hydrated proton in general. 

The protonated water dimer, HsO^, also known as Zundel cation, is the smallest pro- 
tonated water cluster, and has been recently object of intense study. The infrared (IR) 
spectrum of the system has been measured in the gas phase, either using multiphoton disso- 



ciation techniques 0,0] or measuring the vibrational predissociation spectrum of HsO^-RGn 



clusters with RG=Ar,Ne l5|, l8ll. The obtained spectra could not be consistently assigned 
in terms of fundamental frequencies and overtones of harmonic vibrational modes due to 
large-amplitude anharmonic displacements and couplings of the cluster. Hence, more so- 
phisticated theoretical approaches are required. Several theoretical studies have been con- 



ducted over the 
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ast years in order to understand and assign the IR spectrum of the cation 
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The first measurement [6j of the IR multiphoton dissociation spectrum (IRMPD) of H5O2 
spanned the range between 620 and 1900 cm~^. Three main absorptions were discussed 
and assigned, based on a previous quantum-dynamical simulation of the IR absorption 
spectrum on a 4D model of the hydrogen-bond (0-H-O) fragment [q]. Those assignments 
were revisited in the context of newer IRMPD experiments and calculations, producing 
somewhat disparate results , Q, Q- Recent measurements of the IR predissociation 



spectrum of the HsO^ cation in argon-solvate 



and neon- and argon-solvate |5| conditions 



present spectra with a simpler structure than the multiphoton IRMPD ones. It is expected 
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that the spectrum of the HsO^-Nei complex is close to the linear absorption spectrum 
of the bare cation 5|. This spectrum features a doublet structure in the region of 1000 
cm~^ made of two well-defined absorptions at 928 cm~^ and 1047 cm~^. This doublet 
structure was not fully understood until recently, although the highest-energy component 
had been already assigned by Bowman and collaborators to the asymmetric proton-stretch 
fundamental ([0-H-0||]) j^. A similar argument was made by Sauer and collaborators based 
on classical-trajectories calculations |l^. In this respect, it is known from recent classical- 
dynamics simulations on accurate potential energy surfaces (PES) that the [0-H-0||] motion 
features large amplitude displacements strongly coupled to other modes of the system. The 
central-proton displacement would then be involved in most of the lines appearing in the 
IR spectrum, since this motion contributes the largest changes in the dipole moment of the 
13 1 . Recent work by us assigns unambiguously the doublet for the first time fl^ . 



cation 
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which is shown to arise from the coupling of low-frequency motions of the water molecules 
to the central proton motion. 



In this work we undertake the simulation of the IR spectrum anc 



using the multiconfiguration time-dependent Hartree (MCTDH) 16 



dynamics of B.5O2 



17 



18(1 method. In 



doing so we do not rely on a low- dimensional model of the system, but we treat it in its 
full (15D) dimensionality. The obtained IR spectrum is compared to recent experimental 
and theoretical results 5|. The experimental IR spectrum and our results agree well in the 
position and relative intensities of the main spectral features, which allows us to extract 
meaningful conclusions regarding the dynamics of the cluster. All fundamentals and several 
overtones of the vibrational motion of the cation are computed and their properties analyzed. 
The computation of vibrational states in terms of their fully correlated wavefunctions is 
invaluable in the assignment of spectral lines arising from the time-dependent computation 
of the IR spectrum. 

The reported simulations are performed using curvilinear coordinates, which has been 
found crucial for the description of the large amplitude motions and anharmonicities featured 
by the cluster [isl. A description of the derivation of the exact kinetic energy operator (KEO) 
in this set of coordinates and of the representation of the potential energy surface (PES) are 
discussed in the companion paper 15|], hereafter referred to as Paper I. The quality of the 
PES expansion, and the properties of the ground- vibrational state of the system are also 
analyzed in Paper I. The reference PES and dipole moment surfaces (DMS) used are those 
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of Huang et al. ll|, which constitute the most accurate ab initio surfaces available to date 
for this system. 



II. THEORY AND METHODS 



The quantum-dynamical problem is solved using the multiconfiguration time-dependent 



Hartree (MCTDH) method 



11, 



18| . For a brief description of MCTDH see Paper 1. All 



the reported simulations were performed with the Heidelberg MCTDH package of programs 

The definition of the system of curvilinear coordinates used in the description of HsO^ 
is given in Paper I. The 15 internal coordinates are briefly reintroduced here to make the 
discussion self-contained. These are: the distance between the centers of mass of both 
water molecules (R), the position of the central proton with respect to the center of mass of 
the water dimer {x,y,z), the Euler angles defining the relative orientation between the two 
water molecules (waggings: 7a,7_b; rockings: i3a,I3b', internal relative rotation: a) and the 
Jacobi coordinates which account for the particular configuration of each water molecule 
{Ri{A,B), R2{A,B), G{A,B))) where Rix is the distance between the oxygen atom and the center 
of mass of the corresponding H2 fragment, is the H-H distance and 9^ is the angle 
between these two vectors. The grouping of the coordinates into modes is discussed in 
Paper I. Here we repeat only the mode-combination scheme used. The coordinates have 
been grouped into the following five combined modes: Qi = [z,a,x,y], Q2 = [7^,7^], 
Q3 = [R,Uf3^,Ufjg], Q4 = [RiA, R2A, ue^^] and Q5 = [Rib, R2B,U0^g]. 

A. IR Spectrum calculation in the MCTDH framework 



The IR absorption cross-section is given by 20[ |: 



I{E) = :^ z E I (^-l*M,o)P S{E + Eo- E^), (1) 

n 

where I^E^/^.o) is the dipole-operated initial state, i.e., 1"^^^) = ft I'^o), and Eq is the ground- 
state energy. The IR spectrum may be equivalently computed in the time-dependent picture. 
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From Eq. ([T]) follows 

1(E) " 
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Re / e^(^+^°)*/%(t) dt, (2c) 



where the autocorrelation function of the dipole-operated initial state, a^lt), has been de- 
fined in Eq. ([2t) implicitly. The real part, Re, appears because a^{—t) = a*^(t) holds. 
Eq. (jJt) tells that the IR spectrum is obtained in the time-dependent picture by Fourier 
transformation (FT) of the autocorrelation of the dipole-operated initial state. Whether a 
time-dependent or a time-independent approach is more efficient depends on the problem 
at hand. A time-independent approach relies on the accurate computation of the eigen- 
states of the Hamiltonian, {"^n), which may be accomplished by iterative-diagonalization 
methods of the Hamiltonian matrix expressed in some basis, obtaining all the eigenstates 
up to some desired energy, and then using Eq. ([1]) to compute the IR spectrum. The 
method to obtain eigenstates of the Hamiltonian within the MCTDH approach is called 



improved relaxation [2l| and consists essentially on a time-independent multiconfiguration 
self consistent field (MCSCF) computation of an eigenstate, in which both the expansion 
coefficients of the wavefunction and the basis of SPFs (see Eq. (1) in Paper I) are optimized 
to self-consistency. The optimal basis of SPFs varies from state to state, which difficults the 
application of iterative diagonalization procedures to the computation of the whole spec- 
trum. To obtain each excited state with improved relaxation, a new computation is needed 
in which the basis is optimized for that specific state. A comprehensive calculation of an 
IR spectrum by improved relaxation has been recently achieved for 6D problems, namely 



HaCS [2l|, and HONO |22|. The lowest 184 states of A' symmetry of the latter molecule 
were calculated. The improved relaxation approach, however, is not practicable with the 
current computational capabilities to obtain the whole spectrum for a molecule of the size 
of B.5O2 ■ The iterative process of improved relaxation converges only if the space spanned 
by the configurations allows to separate the desired state from neighboring ones. When the 
density of states increases, the number of SPFs necessary to ensure convergence increases as 
well. This makes it impossible to compute states higher in energy. For the present problem 
we could compute states up to about 600 cm^^. However, not only the excitation energy 



but also state coupling determines the feasibility of improved relaxation. For example, we 
could compute the second excited state of the water-water stretch by improved relaxation, 
despite its excitation energy is 1069 cm~^. 

Propagation of a wavepacket by MCTDH, on the other hand, is always feasible. A small 
set of SPFs makes a propagation less accurate, but not impossible. Therefore, the full IR 
spectrum is efficiently calculated in the time-dependent representation. The resolution at 
which different peaks of the spectrum are resolved is given by the Fourier Transform (FT) 
of the damping function g{t) 17|] with which the autocorrelation a^(t) is to be multiplied 
when performing the integral (I2h) to minimize artifacts due to the Gibbs phenomenon. We 
choose g(t) = cos(7rt/2T) and set g{t) = for t > T, where T denotes the length of the 
autocorrelation function. Since we make use of the T/2 trick, a^(t) = {'^*^Q(t /2)\'^ ^Q(t /2)) , 
which holds when the initial wavepacket is real and the Hamiltonian symmetric, T is twice 
the propagation time. The FT of g{t) is known [l^] and its full width at half maximum 
(FWHM) is 

AE = 27, 300 cm~^ fs /T. (3) 

To clearly identify two peaks 100 cm~^ apart, i.e., the doublet structure in the HsOj' spec- 
trum, one needs a resolution of AE = 67 cm~^ and hence a propagation time of T/2 = 200 
fs. Time propagations have been carried out over 500 fs providing a resolution of 30 cm~^. 



Often we make use of the filter diagonalization (FD) technique |23l . |2J, |25|, l26|] to analyze 
the autocorrelation function of the propagated wavepackets, which yields eigenenergies and 
spectral intensities. FD is able to provide accurate eigenenergies from shorter time propaga- 
tions than the ones needed to resolve the IR spectrum. However, FD calculations are more 
sensitive to errors in the autocorrelation function compared to the FT method. Therefore 
these calculations are targeted to obtain accurate energies of a single state or a small group 
of related states from a propagation of an appropriately prepared initial wavepacket. 

The dipole operator is a vector and the operator fi appearing in Eqs. ( fT]l2l) is to be 
interpreted as the scalar product e- fi where e* denotes the polarization vector of the absorbed 
light. As the molecule is assumed to be randomly oriented we average over all orientations 
and obtain 

IiE) = ^{h{E) + 2UE)), (4) 
where Iz and Ix denote the intensities obtained with the x— and 2;— component of the dipole 
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operator, respectively. The factor 2 appears because a^^ = a^y holds. The intensity ly is 
hence not explicitly computed. 

B. Eigenstates of the system and spectrum assignment 

Even if the computation of the IR spectrum in the time-dependent representation is 
feasible, a means of assigning the different lines to specific motions of the system is still 
required. 



Wavefunctions of excited states converged by improved relaxation [2l| contain all the 
possible information of that specific state. The intensity of a given excited state is 
readily obtained by computing the dipole moment | (\l/„,|/i|\l/o) p. But even if an excited 
state of interest has been obtained it is difficult to directly inspect these mathematical 
objects due to their high dimensionality. Moreover, for the higher excited states we do not 
have Y^n) at our disposal but only an autocorrelation function providing spectral lines. In 
both cases we characterize the eigenstates by their overlaps with carefully chosen test states, 
i.e., by the numbers |($test|^n)P- The following procedures are used: 

1. Test states \^test) are generated 

a) by applying some operator O to a previously converged eigenfunction, 

\^test)=Nd\^n), (5) 

where is a normalization constant, e.g., Nz\^q) generates a test state which 
in essence differs from the ground state |\l/o) by a one quantum excitation in the 
proton-transfer coordinate z. 

b) by forming Hartree products, where the SPFs are obtained through diagonal- 
ization of mode-Hamiltonians hj. The hj are low-dimensional Hamiltonians and 
each hj operates on the space of a group of coordinates (see Paper I for details 
on coordinates grouping into modes). Rather than using single Hartree prod- 
ucts one may use linear combinations of products in order to satisfy a symmetry 
constraint. 

2. The overlaps |($iest|^n)P are then computed by 
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a) by direct evaluation of the scalar product if 1^^^) is available. 

b) by Fourier transform of the autocorrelation function 
a(t) = i^tEstl expf— i H t)\^tp.st)- The overlap is obtained via the formula 
3: 

\{^test\^n)\' = ^Re e^^Mt) cos{—)dt (6) 

c) by Fourier transformation of the cross-correlation function 
c(t) = ($fest| exp(— z t) |\E'^^o)- The absolute square of the FT of c(t) at 
energy En must then be divided by the spectral intensities of 1"^^^) to obtain 
|($test|^n)P- This is the fastest method because it does not require additional 
propagations. However, it is also the least accurate procedure. We used this 
method several times to obtain a quick overview, but it was not used to generate 
data reported in this article. 

Excitations related to test states will be denoted {rii qi, n2q2 ■ ■ ■), where ni represents the 
quanta of excitation on coordinate qi (in a separable limit). A test state will be denoted 
by \^niqi,n2q2-.) "whcre the terms with rij = are omitted. In case a test state has been 
generated such that it cannot be represented by this simple notation, it will be defined more 
explicitly. 



III. RESULTS AND DISCUSSION 

FIGURE m AROUND HERE 
FIGURE m AROUND HERE 

The MCTDH spectrum in the full range 0-4000 cm~^ is depicted in Fig. [1] and compared 
to experiment in Fig. [2l The dipole-moment operated ground state /i|\E'o) was propagated 
for 500 fs, yielding an autocorrelation of 1000 fs. The spectrum was calculated according to 
Eq. ([2]) and the FWHM resolution of the spectrum is, according to Eq. ([3]), about 30 cm~^. 

In the following sections the different parts of the spectrum and the peak assignments are 
discussed. TableJT] collects the energies of several states of HsO^ obtained with MCTDH, 
MULTIMODE jsj] and experimental values on the HsO^-Ne system jsl. 

TABLE H AROUND HERE 
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The symmetry labels of the discussed vibrational states are given within the Qie symmetry 
group, which arises due to the feasibility (in the sense of Longet-Higgins [27|) of the wagging 
and internal rotation motions (see Paper I). A character table for the Qiq symmetry group 



is found in Table A-25 in Ref. 



281]. If the internal rotation motion around a is assumed 



to be unfeasible, the symmetry analysis can be performed using the point group. The 
permutation-inversion group Qiq contains the point group as a subgroup, but allows 



additionally to permute the two hydrogens of one of the water monomers [28|, |29|] . The Qiq 
labels reduce to the ones by ignoring (— /+) signs in the labeling provided by Qiq. The 
use of QiQ becomes important when labeling the states related to the internal rotational 
motion. The z component of the dipole moment excites vibrational states of 5^ symmetry, 
while the perpendicular components those of symmetry. 



A. Low-energy region 

The low energy region of the spectrum, below 900 cm~^, has not yet been accessed 
experimentally. The reported spectrum in Fig. [1] shows a strong absorption around 100 
cm~^ followed by two lines of fewer intensity at about 250 and 500 cm~^, respectively. 
The three absorptions arise from excitation of the perpendicular component of the field. 
Modes oscillating at these low frequencies are strongly anharmonic: already in the ground 
vibrational state, the system interconverts between equivalent minima through low barriers 
along the wagging (7^1, 7^) and torsional (a) coordinates 1^. 

FIGURE M AROUND HERE 

The first three excited states associated to the internal rotation (la, 2a, 3a) have been 
computed with improved relaxation and their probability-density along the internal rotation 
a is depicted in Fig. [31 The energy of the (4q;) state has been computed by FD. They 
have excitation energies of 1, 103, 126 and 210 cm~^, and symmetries A]^ , , and A^, 
respectively. Here (la) is the splitting state related to the torsional barrier along a. The 
excitation energy of 1 cm"^ is hence the tunneling split. In fact, the symmetry label of the 
rotationally-split state of a (+) state is obtained by multiplication with A^. 

FIGURE H AROUND HERE 
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FIGURE M AROUND HERE 



Fig. m depicts the probability-density projection on the wagging coordinates for the ground 
vibrational state, which is of symmetry, as well as for one of the two fundamental states 
(wiajWih) of the wagging modes. The fundamental wagging modes are centered at 106 cm~^ 
and belong to the E~ irreducible representation of the Qiq group. Therefore, they are 
doubly degenerate, dark states. The band at about 100 cm~^ (see Fig. [1] center) due to 
the perpendicular component of the dipole corresponds to the states centered at 108 cm~^, 
which belong to the irreducible representation. Such states are combined states of the 
fundamental wagging motion and the (la) state of the internal rotation, of symmetry. 
Thus, both pairs of degenerate wagging states, centered at 106 and 108 cm~^ are split 
states with respect to the internal rotational barrier around a. In a similar way the band 
at about 250 cm~^ corresponds to the combination of the fundamental wagging motion 
and the (3a) internal-rotation state, which belongs to the irreducible representation. 
This combination results in two degenerate states of symmetry. The energies of the 
next three wagging-mode states (^2,^3,^4) are, respectively, 232, 374 and 422 cm~^ and 
they are shown in Figs. \&^, \5jp and Et, respectively. These three states correspond to 
two quanta of excitation in the wagging motions and they can be represented by kets |11), 
(1 20) — |02))/-\/2 and (|20) + |02))/v^, respectively, where the \ab) notation signifies the 
quanta of excitation in the wagging motions of monomer A and B. These states have 
symmetries Bf , B2 and Af , respectively. In the harmonic limit these three states would be 
degenerate. The next 2 states, u'4a,W4b, which are not computed, are degenerate again and 
correspond to kets |31) and |13). State W2 has an energy that nearly doubles the energy of 
the Wix states, since it roughly corresponds to one quantum in state Wia and one quantum 
in state Wib. The strong anharmonicity of the wagging motions as well as the coupling 
between right and left wagging can be further appreciated in the progression of W2, and 
W4 vibrational-state energies. We emphasize again that these three states are degenerate in 
the harmonic limit. In addition, the harmonic- analysis energies of the two lowest wagging- 
fundamentals Wia and Wih are around 300 cm~^ larger than the MCTDH result and do 
not account for their degeneracy, since harmonic normal-modes are constructed taking as 
a reference the C2 absolute minimum. However, as discussed above and in Paper I, HsO^ 
interconverts through low potential barriers between several equivalent minima and has Qiq 
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symmetry. The state ws has four probabihty-density maxima along the 2D space spanned 
by 7a and 7^- They correspond to geometries in which one of the water molecules adopts 
a trigonal-planar geometry (HsO"^ character) and the other adopts a pyramidal geometry 
(H2O character). This state transforms according to the B2 symmetry representation, which 
is also the symmetry of the proton-transfer fundamental. State will play a major role 
due to its strong coupling to the proton-transfer mode. 

The first two fundamentals of the symmetric stretch ([O-Oy], R coordinate) have energies 
of 550 and 1069 cm~^ respectively, and have symmetry, while the rocking fundamentals, 
which are degenerate states, have an energy of 481 cm~^ and are responsible of the band 
appearing slightly below 500 cm~^ in the MCTDH spectrum. In contrast to the wagging 
motion, the rocking motion is fairly harmonic and exhibits only a weak coupling between 
left and right rocking. The energies of the two-quanta rocking states r2, and - defined 
in a similar way to W2, and W4 but with a reverse ordering of the zero-order A;et-states - 
are therefore almost degenerate. These energies read 915, 930 and 943 cm~^. 



B. Doublet at 1000 cm^^ 



The doublet centered at 1000 cm~^ is the most characteristic feature of the IR spectrum 
of HsO^. It is depicted in Fig. [2] (top). The highest energy line has been measured to be at 
1047 cm~^ while the low energy component appears at 928 cm~^ jsl. There is accumulated 
evidence in the literature that the absorption of the proton-transfer fundamental occurs in 
the region of 1000 cm~ ay , [ll|, [12, ll3| . Specifically, the band at 1047 cm ^ in Ref. [5| was 



assigned to the first excitation of the central proton motion [5| based on MULTIMODE [30|] 
calculations. This band is the most intense band of the spectrum since the central proton 
motion along the z axis induces a large change in the dipole-moment of the cation. The low- 
energy component has been recently assigned by us 1^. The doublet is seen to arise from 
coupling between the proton-transfer motion [0-H-0||], the low frequency water-wagging 
modes and the water-water stretching [0-0||] motion. In order to obtain a fundamental 
understanding of the low-energy (|\&[;)) and high energy (I^E^^)) components of the doublet, 
test states were constructed by operating with z on the ground state: l^i^) = z\'^o)N, 
where iV is a normalization constant, and by operating with {R — Rq) on the third excited 
wagging state w^: |$i_r,w3) = {R — Ro)\'^w3)N. Note that \^iz) is characterized by one 
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quantum of excitation in the proton-transfer coordinate whereas \^ri,ws) by one quantum 
in [0-0||] and two quanta in the wagging motion. These two test states were propagated 
and their auto- and crosscorrelation functions were used for FD analysis, which yielded an 
energy of 918 cm^^ for |\E'^) and an energy of 1033 cm~^ for |^^). These energies are in good 
accordance to the peaks in Fig. [2] which arise from the propagation of |^';i,o)- The spectral 
intensities were also obtained by FD analysis. The overlaps of the test states to the states 
making the doublet read: = 0.09, \{^iR^^^\<i/^^)\^ = 0.83 and | p = 0.46, 

I ($ii?,«;3 l^d) P = 0.10. One should take into account that these numbers depend on the 
exact definition of the test states, which is not unique. However, they provide a clear 
picture of the nature of the doublet: the low-energy band has the largest contribution from 
the combination of the symmetric stretch and the third excited wagging (see Figs. \5jp and 
|3j:;), whereas the second largest is the proton-transfer motion. For the high-energy band 
the importance of these two contributions is reversed. Thus, the doublet may be regarded 
as a Fermi-resonance between two zero-order states which are characterized by (1-R, w^) 
and (Iz) excitations, respectively. The reason why the third wagging excitation plays an 
important role in the proton-transfer doublet is understood by inspecting Fig. \5jp and El 
The probability density of state has four maxima, each of which corresponds to a planar 
conformation of H20-H"'" (H3O"'' character) for one of the waters, and a bent conformation 
(H2O character) where a lone-pair H2O orbital forms a hydrogen bond with the central 
proton. When the proton oscillates between the two waters, the two conformations exchange 
their characters accordingly. Thus, the asymmetric wagging mode {w^, 374 cm~^) combines 
with the water-water stretch motion (i?, 550 cm~^) to reach an energy close to the natural 
absorption-frequency of the proton transfer, making these motions coupled. The two states 
of the doublet transform according to the B2 irreducible representation of Qiq. 

One last remark regarding the spectral region of the doublet is the small but noticeable 
absorption, which is appreciated both in the experimental and MCTDH spectra between 
both peaks (see Fig. [2]). This low-intensity absorption is due to the r^, rocking centered at 
930 cm~^ and which belongs to the irreducible representation. In an analogous way to 
the wagging state, the rocking state presents four probability-density maxima along 
the rocking coordinates, each of which consists of one water aligned with the central axis 
while the other is in a bent conformation, and therefore will present some degree of coupling 
to the proton-transfer motion. The low absorption of this band, despite its proximity to 
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the natural absorption-frequency of the proton-transfer motion, is quahtatively explained 
by the fact that the rocking motions do not change the hybridization properties of the water 
monomers and are of a lower amplitude in comparison to the wagging motions. The rocking 
motion hence couples more weakly to the proton-transfer motion and induces a smaller 
change of the dipole-moment as compared to the wagging motions. 

FIGURE M AROUND HERE 
C. 1000-2000 cm-i region 

The region between the proton-transfer doublet and the doublet centered at 1800 cm~^ 
features couplings related to the [0-H-0||] and [O-Oy] motions. The MCTDH spectrum 
reported in Fig. [1] presents three main absorptions in this range, located at 1411, 1741 
and 1898 cm"^. We call the eigenstates producing these peaks I^E'mi), and I'^ms), 

respectively, where the m stands for middle spectral-range. The experimental HsOg'-Ne 
spectrum shows two clearly distinguishable bands at similar positions to the 1741 and 1898 
cm~^ absorptions in the MCTDH spectrum, as depicted in Fig.[2l The spectrum of H502'-Ne 
also shows weak but non- negligible absorption in the region immediately above 1400 cm~^ 
(see Fig. 5 in Ref. 0). 

Propagation of test states followed by Fourier analysis of their autocorrelation functions 
as described in Sec. Ill Bl were used to assign these peaks. The following test states were 
generated using eigenfunctions of low dimensional Hamiltonians: |$i^,i_r), \^iz,2r), \^bu)- 
The test state \^bu) consists on the water-bending with ungerade symmetry, i.e. it is char- 
acterized by (|01) — |10))/-\/2 where the two entries indicate the quanta of bending motion 
of monomer A and B, respectively. 

TABLE M AROUND HERE 

The overlaps where |$) is a test state and |\&) is an eigenstate, are given in Tab. 

[m The analysis contemplates also the states of the doublet at 1000 cm~^, since these states 
are coupled to some extent to the ones in the middle spectral-region, e.g., the eigenstate 
|\Ef^) has a squared overlap of 0.10 with the test state \^bu)- AH states in Tab. [TTlare of B2 
symmetry. 
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State |\l/mi), absorbing at 1411 cm~^ in the MCTDH spectrum, has a largest contribution 
from the \^iz,ir) test state. The experimental HsO^-Ne spectrum shows a weak absorption 
in this region with a lower intensity than the peaks at 1763 and 1878 cm~^, the same trend as 
in the MCTDH spectrum. Based on the general good agreement between the experimental 
and MCTDH spectra we propose that this weak absorption is mainly a combined excitation 
{Iz, IR). A band appearing at 1600 cm~^ in the MM/VCI spectrum jsl was assigned to the 
{Iz, IR) transition. However, the experimental spectrum of HsO^-Ne shows no absorption 
at 1600 cm-i [51. 

State \'^m2), which is responsible for the absorption at 1741 cm~^ in the MCTDH spec- 
trum, has the largest contribution from the \^bu) test state. This peak is then mainly related 
to the ungerade water-bending, and has been already assigned in Ref. ^ and a number of 
works. This peak can be assigned already from a standard normal-modes analysis, since 
its main contribution is from an internal motion of the water monomers, and less from the 
relative motions between them. However, it must not come as a surprise that the eigen- 
state |\&m2) has a total squared overlap of 0.26 with the test states containing one quanta 
of excitation in the proton transfer coordinate, namely |$iz,i_r) and |$iz,2_r)- In a 

fashion similar to the coupling to the wagging motion, as the proton approaches one 
water molecule the equilibrium value of the H-O-H angle shifts to a larger value because this 
water molecule acquires more H3O"'' character. Conversely, the water molecule at a larger 
distance of the central proton acquires H2O character and the angle H-O-H shifts to lower 
values. 

State I^E^ms) is responsible for the absorption at 1898 cm~^ in the MCTDH spectrum, 
the lowest intensity feature of the doublet centered at about 1800 cm~^. This state has the 
largest overlap with the \^iz,2r) test-state. Comparison in position and intensity of this 
peak to the corresponding one in the experimental HsO^-Ne spectrum suggests that the 
latest is mainly related to the {Iz, 2R) excitation. 

Analysis of the values in Tab. HTl shows that the eigenstates in the region 1400-1900 cm~^ 
are characterized by the asymmetric bending and combinations of the proton-transfer fun- 
damental and water-water-stretch fundamental and first overtone, with important couplings 
between them. Such a coupling was already noted by Bowman and collaborators by analyz- 
ing the CI coefficients of their MM/VCI expansion jsj. However, the exact nature of each 
band could not be disentangled, and it was concluded that this "is indicative of large cou- 
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plings among various zero-order states in this region". We show that, despite such couphngs 
exist and play an important role in shaping the spectrum, the different eigenstates involved 
retain each its particular character (note that the largest numbers in Tab. [Ill appear at the 
diagonal) and can be assigned to well defined transitions. 

D. 2000-4000 cm-i region 

Symmetry analysis of the OH stretchings of the water molecules within the Qi^ group 
predicts four vibrational states with labels Af, B2 and E~^, and the corresponding torsional 
splitting A]^ , B2 and E~ states. The states-labeling simplifies to Ai, B2 and E, respectively, 
when adopting the point group, in which the internal rotation around a is treated as 
unfeasible. The use of the smaller group is reasonable here because all the minus- 
states are dark in the linear IR-spectrum. The Ai state corresponds to the symmetrical, 
gerade stretch, the B2 state is the symmetrical, ungerade stretch and the E states are 
the two asymmetrical stretches of gerade and ungerade type. We recall that the notation 
gerade/ ungerade is used to indicate +/— linear combinations of the motions of the two 
monomers, while symmetric/antisymmetric refers to the OH motions within each monomer. 
A symmetry analysis based on the C2 point group, the symmetry group of the absolute 
minimum, labels the symmetric stretches as A and B states and the asymmetric stretches 
again as A and B states. The harmonic analysis results (see Tab [11) yield the symmetric 
stretches separated by less than 10 cm~^ and appear at a lower energy than the asymmetric 
stretches. The asymmetric stretches, despite not being degenerate, are very close in energy, 
separated by less than 1 cm~^. Exact degeneracy of the asymmetric stretches is regained as 
a consequence of the feasibility of the wagging-motions, leading to V2d symmetry. 

The MCTDH spectrum in the water-stretching region reveals two absorptions. The lowest 
energy one is related exclusively to the z component of the dipole, while the largest energy, 
most intense band arises exclusively from the perpendicular component (seen Fig. [1]). The 
lowest energy absorption is related to the symmetric, ungerade stretching, while the highest 
energy absorption is related to the degenerate asymmetric stretchings. The energy separa- 
tion between symmetric and asymmetric stretchin gs i s known to be of about 80 cm~^ as seen 
from different experiments and computations jsl. Isil. so these peaks can be resolved by a 
propagation of about 260 fs. The splitting between symmetric gerade and ungerade stretch- 
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ings is expected to be less than 10 cm~^ according to harmonic results. Thus, a propagation 
of about 2100 fs would be needed to resolve these peaks, which is unfeasible in a reasonable 
amount of computer time for this system. In order to resolve the symmetric stretchings 
two initial test-states were prepared by diagonalization of mode-operators, the symmetric 
gerade and symmetric ungerade (|$sm)) stretchings. Both test states were propa- 

gated and their auto- and crosscorrelation functions analyzed by means of the FD method. 
FD analysis yields the symmetric, in-phase stretching at 3607 cm~^, and the symmetric, 
out-of-phase counterpart at 3614 cm~^. Propagation of a test state \^sa) of the asymmetric, 
doubly degenerate stretchings yields the peak centered at 3689 cm~^. These results are in 
excellent agreement to the bands observed in the experimental HsO^-Ne spectrum jsl. Note 
that the propagation of \^sa) yields a result in better agreement to experiment than the 
peak in Fig. [2] obtained by propagation of the perpendicular component of the dipole. This 
is due to the fact that the wavefunction obtained from application of the dipole operator 
to the ground state is a coherent superposition of several eigenstates that are propagated 
together. The more eigenstates are to be coherently propagated, the more complex is the 
dynamics of the wavepacket and the less accurate becomes the propagation for each individ- 
ual state. For this reason, the most accurate energies (e.g. values in Table [T]) are obtained 
from either improved relaxation to the desired eigenstate, if possible, or propagation of care- 
fully prepared wavepackets (referred to as test states in this work) which have as much 
overlap as possible with the eigenstate of interest or a group of them. As a final remark, 
the spectrum arising from direct excitation of the z coordinate (propagation of the \^iz) 
test state) yields all the peaks in the range 800-2000 cm~^ with almost the same relative 
intensities than compared to the propagation of fizl'^o). However, the propagation of |$iz) 
gives a completely flat spectrum in the region above 3000 cm~^. This indicates that the 
excitation of the stretching motions must occur locally on each water molecule and does not 
depend on the excitation of the central proton. Moreover, the coupling of the stretching 
motions to the central proton is very weak, as indicated by the splitting of only 7 cm~^ 
between gerade and ungerade symmetric stretchings and indicates that the OH stretchings 
of HsO^ are robust to the displacement of the central proton, contrary to the situation seen 
for wagging and bending modes. Despite the excitation of the stretchings is local in each 
water monomer, the calculated band at 3614 cm~^ corresponds to the symmetric, ungerade 
mode of B2 symmetry, since both water molecules are oriented in opposite directions with 
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respect to the ^-direction incident field. The symmetric, gerade mode is a dark state of Ai 
symmetry. 

IV. SUMMARY AND CONCLUSION 

The infrared absorption spectrum of the Zundel cation (HsO^) is calculated by the 
quantum-dynamical multiconfiguration time-dependent Hartree (MCTDH) method in the 
linear-absorption spectral-range - 4000 cm~^. The energies of all fundamentals and several 
overtones related to different motions of the system are reported. The cation is considered in 
its full dimensionality (15D). A curvilinear set of coordinates is used to describe the configu- 
ration of the system. Details on the derivation of the kinetic energy operator used, which is 
exact for total angular momentum J = 0, and on the representation of the potential energy 
surface are given in the companion paper [l^. 

The lowest frequency part of the spectrum, shows a strong absorption at about 100 cm~^ 
due to the combination of the fundamental wagging-modes wi, which are E~ degenerate 
states, and the internal-rotation state (la) of symmetry. The resulting pair of degen- 
erate states is of symmetry and therefore bright. The two absorptions at 250 and 500 
cm^^ are related to the combination state of the fundamental waggings and the (3a) state, 
which results in a pair of degenerate states and to the rocking fundamentals, also of 
symmetry, respectively. These absorptions are related to the component of the field 
perpendicular to the water-water axis. The spectral region at about 1000 cm~^ presents a 
double-peak absorption which is the most characteristic feature of the spectrum. This dou- 
ble peak is seen to arise from the coupling of the proton-transfer motion with a combination 
state involving the wagging mode. The reduced probability density of the wagging 
mode projected onto the wagging coordinates is shown in Fig. \5jp. This state presents four 
probability maxima, each of which corresponds to a water in pyramidal conformation (H2O 
character) while the other is in planar conformation (Hjj'O character), and has an energy 
of 374 cm~^. This state alone absorbs light only very weakly (see Fig. [1]), but the state 
arising from the combination of and the water-water stretching (550 cm~^) reaches an 
energy close to the natural absorption of the proton-transfer at about 1000 cm~^. The 
coupling between these both states explains the doublet absorption at about 1000 cm~^, 
which is interpreted as a Fermi resonance between the combination state (1/2,^3) and the 
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proton-transfer fundamental (1^). 

The region between 1000 and 2000 cm~^ presents three main absorptions at 1411, 1741 
and 1898 cm~^. The peak with the highest intensity in this region is the peak at 1741 cm~^ 
which corresponds to the ungerade bending motion of the water moieties. The eigenstate at 
1741 cm"^ has in addition important contributions from the proton-transfer mode and proton 
transfer combined with the water- water stretching mode (see Tab. [IT]). The eigenstates 
absorbing at 1411 and 1898 cm~^ are both described by one quantum on the proton-transfer 
mode plus one and two quanta excitations, respectively, on the water-water stretching mode. 
The nature of the two eigenstates of the doublet at 1000 cm~^ and the three eigenstates 
between 1000 and 2000 cm~^ is described by clearly defined motions (diagonal elements in 
Tab. ini), e.g. the asymmetric bending. However, they constitute a set of coupled states 
(non-diagonal elements in Tab. [Tll) featuring the wagging, bending, water-water stretch and 
proton-transfer motions. 

The region above 3000 cm~^ presents the direct absorptions of the OH-stretching motions 
starting at about 3600 cm~^. Symmetry analysis of the OH-stretching motions in the or 
QiQ groups reveals that the symmetric gerade and ungerade stretching transform according 
to the Ai and B2 representations, respectively. The ungerade, B2 state absorbs at 3614 cm~^ 
due to the ^-component of the field. The gerade, Ai state, which is dark, has an energy of 
3607 cm~^. The small energy splitting between gerade and ungerade states shows that the 
coupling of these motions to the central proton motion must be very weak. Furthermore, 
our analysis also points out that the excitation of the ungerade, B2 state is completely 
independent of the central proton excitation, and is caused by interaction of the field with 
the local dipole of each monomer. Symmetry analysis reveals also that the asymmetric gerade 
and ungerade states are E degenerate. They absorb at 3689 cm~^ due to the component of 
the field perpendicular to the water-water axis. 



The fact that the HsO^ cation may interconvert between several 



connecting equivalent minima had already been pointed out by Wales 29|], who showed that 



low energy barriers 



the correct symmetry group, because wagging and internal rotation motions are allowed, is 



the permutation- inversion group Qi^ [33(. Our study shows that the symmetry analysis in 
QiQ is necessary to understand some important features of the IR linear absorption spectrum. 
In addition, it has been shown in Paper I that already in the ground vibrational state there is 
non-neglible probability of crossing the internal-rotation barrier. However, the consideration 
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of the cation in the more famihar point group may provide an adequate labelhng of the 
vibrational motions of the system as long as the internal-rotation mode is not involved in 
the considered states. 

The reported calculations are in excellent agreement to the experimental measurements 
of Ref. js] on the predissociation spectrum of HsOj'-Ne. The discrepancy of the MCTDH 
energies reported in Tab. [T] with respect to the position of the measured bright bands lies 
always below 22 cm~^, the average discrepance is 14 cm~^ (see Tab. [T]). Such a remarkable 
consistency between experiment and theory along the whole spectral range represents, on 
the one hand, a validation of the underlying potential eriergy surface of Huang et al. HI] and 
of the mode-based cluster expansion of the potential 15|] used in the quantum-dynamical 
simulations, but is also a clear indication that a suitable set of coordinates was selected to 
tackle the problem [l5|. On the other hand, it provides a validation of the measurements on 
the HsO^-Ne cation by the predissociation technique, since this method gives access to the 
infrared linear-absorption regime without noticeable disturbances caused by the messenger 
atom. The only essential disagreement between experiment and theory is in the intensities. 
The intensities are not measured absolutely but, in comparison to the bright double-peak 
at 1000 cm"^, the structures which appear between 1400 and 2000 cm"^ are too low by 
a factor of 3 when compared to the MCTDH simulation. Such a discrepancy does not 
occur in the spectrum of HsO^-Ar [sl, which displays a relative intensity in that region 
similar to the MCTDH one despite an incorrect shape. Therefore we conclude that more 
investigations are necessary to determine the origin of the discrepancy in relative intensities 
between experiment and theory in some regions of the spectrum. 

The fact that the reported simulations are successful in obtaining accurate results for 
a system of the dimensionality of the protonated water-dimer is to be attributed in great 
part to the MCTDH algorithm, in which not only the expansion coefficients, but also the 
orbitals (here SPFs), are variationally optimal. For a 15-dimensional system the use of only 
4 basis functions per degree of freedom represents of the order of 10^ configurations. The 
largest calculations reported here consist of about 10'' configurations, while already good 
results are obtained by using as few as 10^ configurations (see Tab. 3 in Paper I). Such an 
early convergence of the MCTDH method becomes crucial as high-dimensional problems are 
attempted. 

The reported simulations constitute a new example of the ability of the MCTDH method 
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to tackle high dimensional, complex molecular systems in a rigorous manner, and they open 
exciting perspectives for the simulation and understanding of even more complicated systems 
by means of accurate, non-trivial quantum-dynamical methods. Last but not least, they 
provide important information on the spectroscopy and dynamics of the hydrated proton. 
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TABLE I: Vibrational-excited states as identified in the MCTDH calculations. Comparison is given 
to harmonic-analysis (HO) results in the same surface, the MM/VCI results and the experimental 
results on the HsOg^-Ne cation. In the MCTDH column the subscript D indicates that the state 
was obtained by improved relaxation and the subscript F indicates that the state was identified by 
Fourier analysis. In the |(>I'„|/i|^'o)P column, fi refers to jlz for states and to fix (or fiy) for 
states. In the assignments column, in parenthesis, a number followed by a letter indicates the 
quanta of excitation in that coordinate. Other states are named after their definition in the text. 
As a remainder we note that w, r, b and s refer to wagging, rocking, bending and stretching. When 
meaningful a ket description of the state is given. In the ket description of the OH stretchings, S/A 
indicate symmetric/asymmetric stretching-motion within each monomer, respectively. The last 
column shows the irreducible representation of the Qiq permutation-inversion group to which the 
vibrational state belongs. The irreducible representations of the more familiar point-group 1^2(1, 
which is a subgroup of Qiq, are obtained by simply dropping the upper (-|-/— ) index. 
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TABLE II: Overlaps |(<I>|^)|'^, where \^) are test states and |^') are eigenstates. 
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Figure Captions 

Figure [It Simulated MCTDH spectrum in the range 0-4000 cm~^. Excitation in the z 
direction (top), perpendicular to the 0-H-O axis (middle) and total spectrum, i.e. (l/3)z + 
(2/3)*perpendicular (bottom). Note the different scale of intensities in the perpendicular- 
component plot. Autocorrelation time T = 1000 fs. Absorption is given in absolute scale in 
mega-barns (Mb). (1 Mb = 10"^^ cm^). 

Figure [2t Comparison between the MCTDH spectrum (top) and the HsO^-Ne spectrum of 
Ref. (bottom). The intensity of the experimental spectrum is adjusted in each spectral 
region (800-2000 and 3500-3800 cm~^) using the most intense peak of the MCTDH spectrum 
as a reference. Absorption for the MCTDH spectrum is given in absolute scale in mega-barns 
(Mb). (1 Mb = 10-is cm2). 

Figure |3t Reduced probability density on the internal rotation a for (a) the ground state 
and the first three excited states: (b) la, (c) 2a and (d) 3a. The dotted lines correspond to 
an enlarged scale ( x 10) . The -|- /— symbols are intended to clarify the symmetry properties of 
each state. They indicate the sign of the underlying wavefunction based on a ID computation 
for coordinate a and do not refer directly to the multidimensional wavefunctions from which 
densities are given. 

Figure |4j Reduced probability density on the wagging coordinates ■ja and 7b of (a) the 
ground vibrational state and (b) the first-excited {wia) wagging-mode states. 

Figure [5} Reduced probability density on the wagging coordinates 7^ and 7^ of the excited 
states W2, W3 and W4, characterized by two quanta of excitation. Compare with Table HI 

Figure [6l Schematic representation of the two most important coupled motions responsible 
for the doublet peak at 1000 cm~^ 
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FIG. 1: Vendrell et. al., Journal of Chemical Physics 
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FIG. 2: Vendrell et. al., Journal of Chemical Physics 
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FIG. 3: Vendrell et. al., Journal of Chemical Physics 
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FIG. 4: Vendrell et. al., Journal of Chemical Physics 
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FIG. 6: Vendrell et. al., Journal of Chemical Physics 
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